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SUMMARY 


The  six  degrees  of  freedom  dynamic  equations  of 
aircraft  motion  have  been  documented  in  this  Memorandum  for 
use  in  aircraft  simulations  at  ARL.  Earth  axes  are  chosen 
for  the  integration  of  the  force  quations,  and  body  axes  for 
the  integration  of  the  moment  equations.  The  use  of  quaternions 
to  calculate  aircraft  altitude  and  associated  direction  cosines 
is  described.  This  Memorandum  also  contains  a  brief  description 
of  an  atomspheric  data  subroutine  for  use  in  aircraft  simulation. 
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INTRODUCTION 


In  August  1981,  a  research  program  was  initiated  at  ARL, 
under  Royal  Australian  Navy  sponsorship,  to  study  the  dynamic 
behaviour  of  V/STOL  aircraft.  The  study  will  be  based  on  a  six 
degrees  of  freedom  simulation  of  a  military  jet  V/STOL  aircraft, 
programmed  using  the  Advanced  Continuous  Simulation  Language  (ACSL)* 
on  the  ARL  DEC  System- 10  computer . 

Important  aspects  of  the  simulation  are  the  derivation  of 
the  basic  dynamic  equations  of  aircraft  motion,  and  the  selection 
of  suitable  reference  axes  systems.  Based  on  current  practice, 
earth  axes  were  chosen  for  the  integration  of  the  force  equations, 
and  body  axes  for  the  integration  of  the  moment  equations.  In 
addition,  the  quaternion  parameters  were  chosen  for  the  calculation 
of  the  Euler  angles. 

This  Memorandum  provides  a  documentation  of  the  basic  six 
degrees  of  freedom  dynamic  equations  of  motion  for  ise  in  aircraft 
simulations.  The  axes  systems  used  in  the  set  of  equations  are 
defined  in  section  2,  while  section  3  discusses  the  selection  of 
axes  systems  for  integration.  The  use  of  quaternions  in  determining 
aircraft  attitude  and  performing  axes  transformations  is  dealt  with 
in  sections  4  and  5  respectively.  The  equations  of  motion  are 
derived  and  listed  in  section  6,  and  summarized  in  figure  2.  Section 
7  provides  a  brief  note  on  the  source  of  atmospheric  data  and  lists 
the  atmospheric  models  which  have  been  incorporated. 

2 .  DEFINITION  OF  AXES  SYSTEMS  USED  IN  THE  SET  OF  EQUATIONS 

All  axes  systems  are  assumed  to  be  orthogonal,  right- 
handed  triads,  and  are  shown  in  figure  1. 

2.1  Earth  Axes  (X„,Y„,Z.J 

Ci  £  £ 

The  origin  is  at  a  point  fixed  on  the  earth's  surface, 
typically  at  the  runway  threshold  and  on  the  centreline  (R  <j[_) . 

The  X-axis  points  North;  the  Y-axis  East;  and  the  z**axis  'down' 
toward  the  centre  of  the  earth.  It  is  assumed  that  the  earth  is 
flat  and  non-rotating,  such  that  the  earth  axes  are  regarded  as  an 
inertial  frame. 
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2.2  Body  Axes  (*^»YB'ZB) 

Body  axes  are  fixed  on  the  aircraft  with  the  origin  located 
at  the  aircraft  centre  of  gravity.  The  aircraft  is  assumed  to  be 
rigid,  with  the  X-axis  parallel  to  the  horizontal  fuselage  reference 
line  and  pointing  ' forward ' ;  the  Y-axis  pointing  to  starboard  (right) ; 
and  the  Z-axis  ' downward ' . 

2.3  Stability  Axes  (X_,Y_,Z_) 

bSb 

Stability  axes  are  a  special  set  of  body  axes  used  primarily 
in  the  study  of  small  disturbances  from  a  steady  reference  flight 
condition.  Aerodynamic  data  are  frequently  presented  in  stability 
axes.  These  axes  are  displaced  from  the  body  frame  by  the  angle  of 
attack,  a,  such  that  the  X-axis  in  the  steady-state  is  aligned  with 
the  projection  of  the  relative  v/ind  vector  on  the  aircraft  plane 
of  symmetry;  the  Y-axis  points  to  starboard;  and  the  Z-axis 
•downward* . 


2.4  Air-Path  Axes  (X..,Y,,rZ,J 

"  V*  w  .» 

Air -path  axes  differ  from  body  axes  by  the  angle  of  attack, 
a,  and  the  angle  of  sideslip,  8.  The  transformation  from  body  to 
air-path  axes  is  accomplished  as  shown  in  figure  1,  by  first  pitching 
through  -a,  to  coincide  with  the  stability  axes,  and  then  yawing 
through  8.  The  origin  is  located  at  the  aircraft  centre  of  gravity 
and  the  X-axis  is  aligned  with  the  relative  wind  vector. 

Note :  When  the  wind  velocity  components  are  zero,  the  air-path 

axes  coincide  with  the  flight-path  axes  as  defined  by  Fogarty  and 
Howe  (31 .  Etkin  [2]  refers  to  the  air-path  axes  as  the  air-trajectory 
reference  frame  (or  wind  axes) . 

3 •  SELECTION  OF  AXES  FOR  INTEGRATION  OF  EQUATIONS 

There  are  good  reasons  to  use  mixed  axes  systems  when 
studying  the  many  measurements  and  observations  of  vehicle  motion. 

The  various  options  are  discussed  by  Etkin  [2j  and  other  authors, 
such  as  Fogarty  and  Howe  [31 ,  and  McFarland  [5J . 

3 . 1  Force  Equations 

In  recent  years  there  has  been  a  trend  away  from  the  use 
of  body  axes  for  cue  computation  and  integration  of  force  equations. 
Flight-path  axes  do  have  definite  benefits  in  terms  of  computer 
scaling  and  computation  speed  (3).  However,  where  these  factors 
are  not  of  over-riding  importance,  the  selection  of  earth  axes  is 
becoming  quite  common. 


k 


For  the  proposed  V/STOL  aircraft  studies,  the  earth  axes 
system  was  selected  for  the  computation  and  integration  of  the 
translational  equations  of  motion  for  the  following  reasons: 

(i)  it  gives  commonality  with  the  axes  system  chosen  by  both 
RAE  Bedford  and  NASA  Ames  in  their  simulation  programs 
[1,5]; 

(ii)  it  is  the  logical  choice  for  modelling  an  aircraft  operating 
from  a  moving  ship,  since  there  is  a  need  to  refer  the 
forces  and  motion  of  both  vehicles  to  a  common  reference 
frame ;  and 

(iii)  it  avoids  the  complications  which  exists  with  other  axes 
systems  when  considering  the  presence  of  winds  and 
turbulence,  as  highlighted  by  McFarland  [5], 

3 . 2  Moment  Equations 

The  body  axes  system  is  the  natural  choice  for  the  solution 
of  the  rotational  equations  of  motion  because  of  the  important 
advantage  of  constant  moments  of  inertia  when  calculating  the 
moments  and  angular  motion  of  the  aircraft. 

4.  AIRCRAFT  ATTITUDE  DETERMINATION 

The  attitude  of  an  aircraft  is  defined  in  terms  of  the 
traditional  Euler  angles,  (heading  angle),  0  (pitch  attitude), 
and  $  (roll,  or  'bank'  angle).  In  order  to  avoid  the  problems 
associated  with  the  singularity  in  the  Euler  'rate'  equations, 
which  occurs  when  0=±9O°,  quaternion  components  [7]  or  direction 
cosines  may  be  used  in  the  integration  step. 

Quaternion  components  were  chosen  for  the  following 

reasons : 

(i)  their  time  derivatives  are  always  finite  and  continuous, 
whereas  those  of  the  Euler  angles  possess  singularities ; 

(ii)  the  computations  remain  accurate  as  0  approaches  90°; 

(iii)  it  is  a  four  parameter  method  consisting  of  four  integrations 
with  one  constraint  equation,  whereas  direction  cosines, 
in  principle,  require  nine  integrations  and  six  constraint 
equations. 

The  quaternion  components  are  expressible  in  terms  of  Euler 
angles  as  follows: 
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=  cos4>/2  cos0/2  cosiJj/2  +  sin$/2  sin6/2  sini|>/2 

T  *  sin$/2  cos0/2  cosip/2  -  cos<$/2  sin0/2  sinij*/2 

(1) 

t2  *  cos<$>/2  sin0/2  cosiJ//2  +  sin$/2  cos0/2  sinip/2 

■  cos$/2  cos0/2  sin^/2  -  sin<ji/2  sin0/2  cosip/2 

The  quaternion  component  time  derivatives  are  given  by, 

tq  «  -  1/2  (Pt1  +  Qt2  +  Rt3) 

t  »  1/2  <Pt  -  Qx  +  Rt.) 

*  (2) 

t2  *  1/2  (Pt3  +  Qtq  -  Rt1) 

t3  -  -  1/2  (Pt2  -  0TX  -  Rtq) 

where  P,Q,R  are  angular  velocity  components  about  body  axes,  and 
V  +  T1*  +  T2a  +  T3a  =  1  (3) 


Failure  to  normalize  the  quaternion  components  at  each  iteration  can 
result  ir  the  integration  becoming  unstable. 


Euler  angles  may  be  derived  from  the  quaternion  components 
by  using  the  following  relationships, 


$  =  tan 


T2T3  +  Vl 

V+V"1'2 


=  tan 


•1  ’V 

JSJ 


9  =  tan 


0  *  tan 


■1  _ T0T2  ~  T1T3 _ 1  .  -1  ~L3 

f  <  v+V"1/2)  s+<t1t2+t0t3)  *ri7arj  * tan  [TIp+EpP72 


■1  I  T1T2  *  T0T3  | 

LV+V‘1/2J 


-i  rv 

_li. 


The  initial  Euler  angles  are  used  to  determine  the  initial 
quaternion  components,  which  are  in  turn  used  to  calculate  the  direction 
cosine  parameters  for  use  in  axes  transformation  computations.  The 
quaternion  components  are  updated  at  each  iteration,  using  equation 
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(2) ,  such  that  the  direction  cosines  are  recalculated  for  use  in 
the  equations  of  motion,  while  the  Euler  angles  are  calculated  as 
cutput  data  only. 


5.  AXES  TRANSFORMATION 


Transformation  of  a  set  of  variables  from  body  axes  to 
earth  axes  (or  vice-versa)  is  conveniently  achieved  by  use  of  direction 
cosines  [63 ,  which  are  obtained  in  terms  of  the  quaternion  components 
by  the  following  relationships, 


h 

L 

L 


1 

2 

3 


2  (V  +  Tia)  • 

2  (T1T2  +  T0T3> 


2  (T1T3  '  T0T2) 


M1  =  2  (T1T2  “  T0T3> 

m2  =  2  <V  +  V>  '  1 


N 

N 

N 


1 

2 

3 


2  (V1T3  t  W 
2  (T2T3  '  Vl) 

2  (T0a  +  T32)  "  1 


(7) 


6.  EQUATIONS  OF  MOTION 

Figure  2  is  a  summary  of  the  overall  six  degrees  of  freedom 
dynamic  equations  for  the  case  of  a  flat  earth. 

6 . 1  Force  Equations 

The  aerodynamic  force  components  are  frequently  computed 
along  stability  axes,  and  the  propulsive  force  components  are  usually 
supplied  in  body  axes.  The  resolutes  along  stability  axes  are: 

=  X„  cosa  +  Z„  sina  +  X 

P  P  A 

"  YP  +  *A  (8) 

**  Z  cosa  -  X_  sina  +  Z„ 

P  P  A 


XS 


YS 


zs 


i 
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The  total  force  components  in  earth  axes  are  obtained  by 
the  transformation  of  the  forces  in  stability  axes,  using  direction 
cosines,  and  the  addition  of  the  gravitational  force  component,  as 
shown  by  equation  {9).  ’g’  is  assumed  constant  such  that  the 
calculated  altitude  in  figure  2  is  the  geopotential  height,  as  used 
in  standard  atmosphere  calculations. 


?XE  * 

(I^cosa 

+ 

II,  sina) 

Fxs 

+ 

-Vvs 

+ 

(Mucosa 

-  I^sina) 

F 

ZS 

fye  * 

(L^cosa 

+ 

D^sina) 

Fxs 

+ 

M2FYS 

+ 

(N^coso 

-  L2sina) 

Fzs  (9) 

F2E  " 

(L^cosa 

+ 

N^sina) 

Fxs 

+ 

M3FYS 

+ 

(N^cosa 

-  L3sino) 

F_„  +  VJ 

ZS 

The  components 
are  obtained  simply  from 

of  acceleration  with  respect  to  the  earth 
the  dynamic  equations: 

WE 

= 

FXE/m 

^EE 

= 

'V" 

(10) 

*DE 

= 

where  ir  is  the  aircraft  mass,  and  virtual  mass  effects  are  ignored. 


The  velocity  components  of  the  aircraft  relative  to  the 
earth  are  obtained  from  the  direct  integration  of  equation  (10). 

The  components  of  the  wind  velocity  relative  to  the  ground  are  then 
added,  giving  the  velocity  components  of  the  aircraft  relative  to 


(11) 


►  f 

V 

-  V 

NE 

WN 

V 

-  V 

EE 

WE 

V 

-  V 

DE 

WD 

The  components  of  the  aircraft  velocity  vector  relative  to 
the  air,  in  body  axes,  are  derived  using  equations  (7)  and  (11), 
as  shown  by  equation  (12). 


U6  -  ft,  V  .  ft,  VE)  ♦  (L3  vc> 
VB  -  (Mj  VB>  ♦  (t,2  VE)  ,  (M3  Vo) 


«B  -  (Nj  V  +  (N2  Ve)  ,  (»3  V[)l 


(12) 


-7- 


The  resultant  velocity  (or  True  Airspeed)  of  the  aircraft 


is  then: 


V  =  (U  *  +  V  2  +  w 

TAS  B  B  B  ' 


The  angles  of  attack,  a,  and  sideslip,  8,  are  also 
determined  using  equations  (12)  to  (15); 


o  = 

tan  1  | 

[v°b1‘ 

in  the  range 

+1806 

(14) 

8  = 

tan  *  j 

:  w 

*  v>1/2] 1 

in  the  range  ±90° 

(15) 

8  = 

.  -1 

sin 

[VW 

] 

(16) 

Time  derivatives  of  these  angles  are  reouired  in  the 
aerodynamic  force  and  moment  calculations.  Adequate  estimates  can 
be  obtained  from  polynomial  fits  to  the  current  values  and  those 
from  a  few  preceeding  time  steps.  Tomlinson  [8]  uses  a  quadratic  fit 
to  estimate  a  and  8. 

e.g.  d(t)  -  a(t~2At)  -  4a(t-At)  +  3a (t)  (17) 

2At 

Plight-path  parameters  are  derived  directly  from  the  earth 
axes  velocity  vector  components  using  equations  (18)  to  (20) : 

Ground  speed,  . 

V  =  (v  a  +  V  *) 1/2  (18) 

GM  '  N'E  EE  ' 

Flight-path  angle, 

T  -  t“'1  Ke'VI  ,19> 

Angle  of  track,  east  of  north 

X  -  (VV)  (20) 

The  positional  coordinates  of  the  aircraft's  centre  of 
gravity  are  derived  by  integrating  the  earth  axes  velocity  vector 
components  ■, 

Distance  travelled  Worth, 

x  =  V  (21) 


Distance  travelled  East, 
y  =  v 


Altitude , 


Af.T  =  -v 

DE 


(23) 


The  additional  parameter  of  aircraft  range  is  calculated 
using  the  definition: 

Rg  -  <x*  +  y*)1/2  (24) 

Linear  accelerations  at  the  aircraft's  centre  of  gravity 
(in  units  of  'g'),  are  calculated  by  transforming  the  linear 
accelerations  in  earth  axes  (equation  (10))  to  body  axes,  using 
equation  (7),  and  dividinq  by  'g'.  Linear  accelerations  at  the 
pilot's  head  and  elsewhere  (e.g.  for  accelerometers)  can  then  be 
derived  as  required  by  the  user. 

6 . 2  Moment  Equations 

The  total  moments  acting  on  an  aircraft  consist  of  aero 
dynamic  and  powerplant  components.  The  powerplant  components  whi 
include  gyroscopic  moments  due  to  powerplant  rotors,  and  thrust 
alignment  moments,  are  normally  given  in  body  axes.  The  aerodynamic 
moments  are,  like  the  aerodynamic  forces,  frequently  given  in  stability 
axes  [3] .  If  the  aerodynamic  moments  are  given  in  body  axes,  the 
simplification  of  equation  (25)  is  obvious  (i.e.  set  a=0) . 

L  =  L  coso  -  W  sina  +  L„ 

A  A  P 

M  -  Mft  +  Mp  (25) 

N  =  L  sina  +  coso  +  N 

A  A  P 

If  it  is  assumed  that  the  aircraft  has  a  plane  of  symmetry, 
such  that  the  products  of  inertia  I  and  I  arc  zero,  then  the  body 
axes  angular  accelerations  can  be  calculateaby  using  equation  (26) . 

p  »  L.ca+  ::.c2  +  <p.c3  +  r.c^)  q 

Q  *  U.C5  +  <R*-P»)  C,  +  R.P.C?  (26) 

R  =  tJ.C8  +  L.C2  +  (P.C9  -  R.C3)  Q 


whore 
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Ixx  rzz  Ix>’3 

X7Z/C0 

Ixz/co 


C2  (IXX 


+  ^ 


C1  (IYY 


I  J  -  Cl 


Z2 


2  X£ 


=  1/1 


YY 


C5IXZ 


C5  (IZZ 


IXX) 


^O 


ce  (Ixx 


W  +  C2rxz 


The  constants  CQ  to  Cg  are  evaluated  during  initialization. 

The  aircraft  angular  velocity  components  may  then  be 
obtained  by  integrating  equation  (26) . 

7.  ATMOSPHERIC  DATA 

Ambient  atmospheric  temperature  (Tamb) ,  pressure  (Pamb) , 
and  density  (p) ,  are  obtained  as  a  function  of  altitude  from  a 
subroutine  version  of  the  program,  ATMOS. PHH. 

ATMOS. PHH  calculates  the  variation  of  atmospheric  properties 
with  altitude  and  has  an  altitude  range  from  zero  to  100,000  feet. 
There  are  six  options  available  to  the  user;  selection  being  made 
by  setting  the  flag  "KEYAIR"  in  accordance  with  Table  1. 

ATi'OS.PHH  can  be  easily  modified  to  calculate  temperature/ 
pressure /density  ratios,  Mach  number  and  speed  of  sound  if  required. 
The  criteria  required  to  calculate  the  ARDU  Tropical  Atmosphere  are 
those  presented  by  Kipp  (4] . 


k 
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TABLE  1: 

ATMOSPHERE  OPTIONS  OF  ATMOS. PHH  ! 

J 

KEYAIR 

ATMOSPHERE 

DATA  REQUIRED 

1 (default) 

ICAO  Standard  Atmosphere 

ALT 

2 

ICAO  Sea  Level  conditions  at 

- 

all  times 

3 

Off-Standard  ICAO 

ALT , TDAY , QNH , 

HAFR 

4 

ARDU  Tropical  Atmosphere 

ALT 

5 

ARDU  Sea  Level  conditions  at 

- 

all  tires 

6 

Off-Standard  ARDU 

ALT , TDAY , QNH , 

HAFR 

-  i 

8.  CONCLUSION 

The  basic  six  degrees  of  freedom  dynamic  equations  of 
aircraft  motion  have  been  documented  in  this  Memorandum  for  use  in 
aircraft  simulations  at  ARL.  The  earth  axes  system  was  selected  for 
the  integration  of  the  force  equations,  and  the  body  axes  system  for 
the  integration  of  the  moment  equations.  Quaternion  components  have 
been  used  to  calculate  Euler  angles  and  direction  cosines.  A  routine 
for  calculating  atmospheric  variables  has  been  provided. 
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NOTATION 


ALT 

Ci  (i*0,9) 

F  P  P 
XE'  YE'  Z £ 

fxs'fys,fzs 

g 

HAFR 


Vyy'Vz^xz 

KEYAIR 

L,M,N 

Vma'na 

VVKP 

Li,Mi,Ni  ( i— 1 ,3) 

MSL 

m 

Pamb 

P,Q,R 


QNH 

rg 

Rxt 

Tar.ib 

TDAY 


VGR 

VN'VE'VD 


vne,vee'vde 


TAS 

VHH'VMB'VWD 


altitude  (ft,  m) 
inertia  constants 

applied  forces  in  earth  axes  (lbf,  N) 

applied  forces  in  stability  axes  (lbf,  N) 

gravitational  constant  (32.174  ft/s2,  9.81  m/s2) 

height  of  the  airfield  reference  point  above 
MSL  (ft,  m) 

moments  and  product  of  inertia  (lb  ft2.  Kg  m2) 
flag  to  select  desired  atmosphere  in  ATMOS. PHH 
total  moments  (ft-lbf,  N-m) 
aerodynamic  moments  (ft-lbf,  N-m) 
powerplant  moments  (ft-lbf,  N-m) 
direction  cosines 
mean  sea  level 
aircraft  mass  (lb.  Kg) 

ambient  atmospheric  pressure  (lbf /ft2,  N/m2) 

angular  velocity  components  about  body  axes 
(rad/s) 

airfield  pressure  altitude  above  MSL  (millibars) 

aircraft  range  (ft,  m) 

centreline  of  runway  threshold 

ambient  atmospheric  temperature  (C) 

ambient  air  temperature  at  airfield  reference 
point  (C) 

body  axes  velocity  components  (ft/s,  m/s) 

aircraft  ground  speed  (ft/s,  m/s) 

components  of  airspeed  in  earth  axes  (north, 
east,  down)  (ft/s,  m/s) 

components  of  velocity  relative  to  earth 
(ft/s,  m/s) 

true  airspeed  (ft/s,  m/s) 

components  of  total  wind  velocity  relative  to 
earth  (ft/s,  m/s) 

aircraft  weight  i.e.  W  =  mg  (lbf,  N) 
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NOTATION  (CONT.) 


XB,Y  ,Za 
AAA 


VV2b 

VVZE 

W*P 

xs'Ys'2s 

VVzw 

x»y 

a 

6 

Y 

6  ,6  ,6 
a  e  r 


p 

Ti  (i=0, 3) 
X 


aerodynamic  forces  (drag,  sideforce,  lift) 

(lbf ,  r.) 

body  axes  reference  frame 
earth  axes  reference  frame 
propulsive  forces  (lbf,  ”) 
stability  axes  reference  frame 
air-path  axes  reference  frame 
positional  coordinates  (ft,  m) 
angle  of  attack  (rad.) 
angle  of  sideslip  (rad.) 

flight-path  angle,  or  ’angle  of  climb’  (rad.) 

control  surface  deflections:  aileron,  elevator 
rudder  (rad.) 

Euler  angles  of  pitch,  roll  (bank)  and  yaw 
(heading)  (rad.) 

air  density  (lb/ft3,  Kg/m3) 

quaternion  components 

angle  of  track,  east  of  north  (rad.) 


Subscripts 


amb 

A 

B 

E 

P 

S 

K 


VVWD 


ambient 

aerodynamic  contribution 
body  axes 
earth  axes 

powerplant  contribution 
stability  axes 
air-path  axes 

North,  East  and  Downwards  wind  velocity 


A  dot  over  a  variable  denotes  the  first  derivative  with  respect  to 
time. 


FIG.  1:  RELATIONSHIP  BETWEEN  EARTH  AXES.  BODY  AXES,  STABILITY  AXES 

AND  FLIGHT-PATH  AXES 


Propulsion 

Forces  Resolution  of  forces  into  Stability  Axes 

Xp  Fxs  =  Xpcosa  +  Zp  sin  a  +  XA 


Stability 

Axes 

Forces 


Fx  e  —  (L^  cos  fl  ^  N,  sin  a)  F  Xg  +  M^F  yg 


F  y  e  =  (L2cosa  +  N2sina)Fxs  +  M2FVS 


Control  Surface 
Displacements 


Lp  y  Mp  |  NP 
Propulsion  Moments 


Up-dated 


Normalized 


L,  =  2 


Quaternions  Normalize  Quaternions  Quaternions!  1^-2 


L3  -  2 


M,  =  2 


1)FXS  +  M,FYg  +  (M,  cosa  -  L,  sina)Fzs 

)  F  xs  +  M2FYS  +  <N2cosa  -  L2sina)Fzs 

)  Fxs  +  M3FYg  +  <N3  cosa  —  L3  sin  a)  Fzs  +  W 


Vne 

111 

u.* 

II 

/m 

Vee 

UJ 

u> 

a 

/m 

Vde 

=  fze 

/m 

Earth  Axes 
Velocity  Components 


P  Q  R 

_ T  1  i 

.C,  +N.c2+  (PC,  +  R.C4)Q 
I.C5  +  (R;  —  P: )  Cg  +  R.P.Cj 
LC8  +  L.Cj  +  (P.C9  -  R.Cj)  Q 


Body  Axes 
Angular  Rates 


Wind  Velocity 
Components 

I  ^WN  I  VWE  I  VWD 


VN  ~  VNE  -  VWN 

vEE 

VE  =  VEE  _  VWE 

VDE 

VD  =  VDE  ~  VWD 

VN|  VE|  vd| 


=  IVNE.’  +  VEEJ 


-  tan  ( -  Vrp  V 


Tan  [Vet  V*. 


X 

=  VNE 

V 

VFE 

alt 

~VDF 

Directional 
Cosi  nes 

h  =  21  r0:  +  r,5  )  -  1  L13 

Lj  =  2I-,  r- +r0r3) 

M3 

L3  =  2(7^3  -r0r2)  - 

Mi  =  2(r,r2  -  7073!  N3 

M  2  —  2  (  Tq  ■+■  t  2  )  1 

Mj  =  2(72  73  +  707,) 

m,,m2 

N1  =  2(7,  73  +7q72  )  - 

Nj  =  2(7273  -  707,l  N,  ,  N2 
N3  =  2(70,+73M  -  1 
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Velocity  Components 
Relative  to  Air 
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tan 

1  l-L3/(L,-  +  L2:)*l 
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1  [M3/N3] 

c 

CD 

•'  (L1  VN  1 

+  <L2Ve  ) 

*  IL3 Vp  1 

VB 

-  (M,  VN  1  +  (M2Ve 
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WB 
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Body  Axes 

Velocity 

Components 


UB  VB  WB 


Il,.3  I  M,  3  I IM, 


^ 

Euler  Angles 


VTAS  '  <UB'  +  VB‘  +  WB-  I  ■ 


a  -  tan  1  [WB  UBI 


j3  =  tan'1  |VB  (UB:  ♦  WB:I  i 


OCK  DIAGRAM  OF  COMBINED  EARTH  AXES/BODY  AXES  SYSTEM  (FLAT  EARTH) 

FOR  AIRCRAFT  MOTION. 
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